Frequency and surface dependence of the mechanical loss in fused silica 
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We have compiled measurements of the mechanical loss in fused silica from samples spanning a 
wide range of geometries and resonant frequency in order to model the known variation of the loss 
with frequency and surface-to-volume ratio. This improved understanding of the mechanical loss 
has contributed significantly to the design of advanced interferometric gravitational wave detectors, 
which require ultra-low loss materials for their test mass mirrors. 
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INTRODUCTION 

As part of the research and development for the 
LIGO 1] and TAMA j3] gravitational wave detectors, we 
have conducted investigations into the internal friction 
of fused silica. Displacement of the interferometer's mir- 
ror faces arising from thermal motion of the fused silica 
test mass mirrors sets a fundamental limit to the detec- 
tor sensitivity. The frequency distribution of this noise 
is directly related to the internal friction of the mirror 
material. 

An Advanced LIGO detector has recently been pro- 
posed with better sensitivity than initial LIGO. The 
Advanced LIGO mirror thermal noise must be as low 
as possible. Two materials have been under consider- 
ation for the mirror substrate: fused silica and single 
crystal sapphire. To its advantage, sapphire has the 
higher Young's modulus and a low bulk mechanical loss 
(0 < 3 X 10-9) 0. However, sapphire also has high 
thermoelastic noise j^. 

In the advanced detectors, thermal noise in the mir- 
ror coatings makes a significant contribution to the total 
noise budget in the central frequency region of 30-500 Hz. 
Discussion on the mechanical loss in the mirror coatings 
can be found elsewhere 13 ■ 

Recent measurements on the mechanical loss in fused 
silica have revealed a dependence on frequency 0,^3 ^'^'^ 
on surface-to- volume ratio 0, ITsf . This paper com- 
bines data from several of these research groups in order 
to model both of these effects. The frequency depen- 
dence j3£the loss agrees well with results from Weidersich 
et al. 14] . In that work, loss data spanning six decades in 
frequency is modeled by an asymmetric double-well po- 
tential in the bond angle. Together these results provide 
a more complete picture of the loss in ultrapure glasses 
and a more physically motivated prediction for the ther- 



mal noise in advanced interferometric detectors. It was 
previously predicted that fused silica's loss dependence 
would make it suitable for low frequency detectors (10 - 
100 Hz) jl^- Indeed, this model's prediction of a very 
low mechanical loss in the LIGO frequency regime has 
motivated the recent selection of fused silica as the Ad- 
vanced LIGO test mass substrate [l6j. 



THEORY OF LOSS IN FUSED SILICA 

The thermal noise motion of the mirror surface is re- 
lated to the internal friction of the substrate by the 
fluctuation-dissipation theorem The internal fric- 

tion of very pure fused silica is associated with strained 
Si-O-Si bonds, where the energy of the bond has min- 
ima at two different bond angles, forming an asymmetric 
double-well potential. Redistribution of the bond angles 
in response to an applied strain leads to mechanical dis- 
sipation, which at audio frequencies has a peak in the 
cryogenic range 20-60K. Because fused silica is an amor- 
phous material, there is a distribution of potentials which 
must be inferred from measurements of the dissipation. 
It can be shown 0| that the frequency dependence of 
the loss should exhibit a power law spectrum with expo- 
nent ksT /Va at low temperatures. Both this power law, 
with Vb/^B = 319iC, and the distribution of potentials 
have been measured |l4| . The power law exponent of a 
relaxation process cannot exceed 1, and is expected to 
saturate near 300 K. At room temperature the exponent 
is 0.76. 

At elevated temperatures there is another loss peak 
arising from a double-well potential associated with the 
Si-O-Si bond angles. For this peak the bond angle shift 
and potential barrier are much larger; the double-well of 
the cryogenic loss peak is a small feature at the minima 
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of this larger potential well. At room temperature, ther- 
mal fluctuations allow the bonds to span the cryogenic 
double-well but not to cross the larger potential barrier, 
where Vb/fc^ = 3.54 x lO^ifl^. The calculated internal 
friction for this loss peak at audio frequencies and room 
temperature is utterly negligible compared to other loss 
mechanisms cited herein. 

A separate loss mechanism exists in the surface of the 
glass. The contribution from the surface loss depends on 
the mode of the sample. The total energy lost per oscil- 
lation in an isotropic sample undergoing slowly decaying 
vibration, can be described by the integral of the local 
loss angle, (/)(r} with the energy density Psi^) 

AE ^ 2iT p^{r) cj){r) d^'r (1) 

where V is the sample volume. Assuming that the lo- 
cal loss angle is constant and equal to ij^buik everywhere 
except within a distance h of the surface, and that the 
energy density in that surface layer is approximately the 
energy density at the surface, then the loss can be ex- 
pressed as [13 

= 0buik + A* "sf (2) 

where S is the surface area of the sample and /i is a factor 
of order unity that depends on the mode shape. The 
surface loss parameter, as, is typically several picometers 
for flame polished or flame drawn fused silica but much 
higher for abrasively polished surfaces. 

EXPERIMENTAL METHOD 

The measurements at Syracuse University (SU) [lolllll 
[l^ Il3l | were performed on fiber/rod samples with diam- 
eters ranging from 0.1-8 mm over resonant frequencies 
less than 5 kHz. The samples were drawn from and left 
attached to a massive bob of Suprasil 0, thus forming 
a cantilever beam. This bob was welded to a vibration 
isolating suspension formed by similar silica bobs con- 
nected by thin silica fibers. In a vacuum of w 10^^ torr, 
the samples were made resonant by an electrostatic comb 
exciter, and their position was measured using a shadow 
sensor. 

The measurements at University of Tokyo 9] were per- 
formed on cylindrical samples with optically polished sur- 
faces. The diameters and heights were 70 mm and 60 
mm, respectively. The samples were annealed in a vac- 
uum furnace. To exclude the support loss, the samples 
were supported at nodal points of their vibrational modes 
during the Q measurements. 

The Caltech measurements were performed on a spare 
input test mass for the initial LIGO interferometers, a su- 
perpolished right cylinder made from Suprasil 312 with 
a diameter of 25.4 cm and a thickness of 10 cm. It was 



suspended in a « 10~^ Torr vacuum by a loop of pol- 
ished stainless steel wire greased with lard. The elastic 
modes of the mass were excited with an electrostatic ac- 
tuator and the mode amplitude was monitored using a 
birefringence sensor. Since friction at the wire could re- 
duce Q, only modes with small motion at the point of 
wire contact were used in the fit. 



MODELING METHOD 

Resonant Q measurements from each of the labs were 
submitted for generating this model of the loss. The 
measurements spanned several types of fused silica, V/S 
ratios from 0.03 - 28 mm, and frequency up to 10^ Hz. 
The data was first separated by silica type since the loss 
is known to vary significantly between varieties of fused 
silica 0,0, 1131 ■ Only Suprasil 2 and 312 had sufficient 
data to warrant a fit over both frequency and V/S ratio. 
Characteristics of these samples are listed in Table 



Label 


Type 


D 


h 


V/S 


Surface 


Anneal 


Lab 


PI 


312, cyl. 


254 


100 


28 


SP 


None 


Caltech 


K12 


312, cyl. 


70 


60 


11 


SP 


980°C, vac. 


Tokyo 


SU2 


312, rod 


3 




0.75 


FP 


1025°C in Ar 


SU 


SV4 


312, rod 


8 




2 


FP 


950°C in Ar 


SU 


AGS 


2, rod 


3.5 


188 


0.88 


FP 


None 


SU 


AHl 


2, rod 


0.300 


108 


0.075 


FP 


None 


SU 


ANl 


2, rod 


0.318 


160 


0.080 


FP 


None 


SU 


ABl 


2, rod 


0.062 


175 


0.016 


FP 


None 


SU 


ACl 


2, rod 


0.340 


310 


0.085 


FP 


None 


SU 


AFl 


2, rod 


0.120 


130 


0.030 


FP 


None 


SU 


K13 


2, cyl. 


70 


60 


11 


SP 


900°C, vac. 


Tokyo 



TABLE I: Sample characteristics: Type lists Heraeus Suprasil 
variety and shape. Samples are cylinders with diameter (D), 
height (h), and volume-to-surface ratio (V/S) given in mm. 
Surface types are superpolished (SP) and flame polished (FP). 

We chose a model for the mechanical loss that included 
terms describing the frequency dependence, the surface 
loss, and the thermoelastic loss. The loss function took 
the form: 

V 

4'if,-^) = "/"surf + 0bulk + </'th (3) 
= Cli^y + C2{f/l + CiCl^ti, (4) 

where Ci ~ pas from Eqn. |2] Given that the surface 
loss term only contributes significantly to the rod (fiber) 
samples, we have assumed for all samples that p ^ 2 
which is appropriate for cylindrical rods. We have also 
not distinguished the loss angle arising from the Young's 
modulus from that due to the shear modulus. 
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FIG. 1: Suprasil 2 mechanical loss data: Best fit surface (up- 
per) and Deviation in units of sample variance (lower). 



FIG. 2: Suprasil 312 mechanical loss data: Best fit surface 
(upper) and Deviation in units of sample variance (lower). 



Type 


Ci (pm) 


Ca (xlO"") 


C3 


C4 


2 


12.1 ± 0.8 


1.18 ± 0.04 


0.77 ± 0.02 


0.61 ± 0.05 


312 


6.5 ± 0.2 


0.76 ± 0.02 


0.77 ± 0.02 





TABLE II: Fit coefficients for Suprasil 2 and 312. 



The thermoelastic loss, </)th, which is neghgible in all 
but the thinnest fiber samples, is described for fibers by: 



Ya'T 



pC„, l+(27r/r)2 



(5) 



{(P pC^)/ {13.55 k) 



where Y is the Young's modulus, a is the coefficient of 
thermal expansion, T is temperature, p is the density, Cm 
is the mass specific heat capacity, d is the diameter, and 
K is the thermal conductivity. We fit the amplitude of (f>th 
to account for small changes in the coefficient of thermal 
expansion among samples. Variations in the fiber diam- 
eter can also slightly alter the shape of the thermoelastic 
peak. Neither of these effects significantly affect the fre- 
quency or surface loss terms. 



Measurements of large resonant Q's are subject to nu- 
merous mechanisms that can greatly reduce the Q and 
few processes that can increase it. These effects produce 
a distribution in the systematic error that is asymmet- 
ric, heavily skewed toward lower and unique for each 
experiment. Standard data analysis techniques based on 
normally distributed error, such as linear least squares 
(LLS) fitting, are therefore inappropriate for analyzing 
our full data set. We circumvent this problem by first lim- 
iting our data to the best measurement at each (/, V/S) 
point for each sample. A LLS fitting routine is ap- 
plied with the sample variance approximating the ac- 
tual variance of the data. This method is commonly 
used in analyzing mechanical loss measurements where 
the lowest loss measurement closely approximates actual 
mechanical loss for a sufficiently large set of measure- 
ments The results of the method are displayed 
in Figurcnifor Suprasil 2 and in FigureHlfor Suprasil 312. 
The fit coefficients are listed in Table Hll The frequency 
dependence, C3, agrees well with results from Weidersich 
et al. 01 . The thermoelastic amplitude, C4, is similar to 
earlier measurements |loj . Assuming no unforeseen loss 
mechanisms, the Advanced LIGO test masses (V/S « 40 
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mm) have a predicted loss (0(100 Hz) « 4 x 10 that 
is a several-fold improvement over previous estimates. 
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FIG. 3: Estimated Advanced LIGO thermal noise for a 
Suprasil 312 test mass substrate and for two mirror coatings: 
the best measured and the research goal. Laser quantum noise 
provided for comparison. 



Coating Loss 


Loss Angle 0| | 


BNSI Range 


BH/BH Range 


Measured 
Goal 


1.6 X 10"^ 
5.0 X lO"-'' 


190 MPc 
230 MPc 


840 MPc 
1060 MPc 



TABLE III: The distance a single Advanced LIGO interfer- 
ometer could detect a neutron star or 10 Mq black hole binary 
inspiral, assuming a Suprasil 312 test mass and two possible 
mirror coatings: the best measured and the research goal. 



measurements to date |23j and the research goal. The 
sensitivity goal for a single Advanced LIGO interferom- 
eter is to observe BNSI, averaged over sky position and 
polarization, to a distance of « 200 Mpc. (See Harry 
for a description of a LIGO range calculation.) 



CONCLUSIONS 

We have shown that the mechanical loss of fused silica 
can be described by a model that includes surface loss and 
a frequency dependent bulk loss. The frequency depen- 
dent loss, thought to arise from an asymmetric double- 
well potential of the bond angle, agrees well with earlier 
measurements [l^ that spanned six decades in frequency. 
This improved understanding of the loss indicates that 
at large geometries and low frequency, fused silica is an 
excellent material for test masses in advanced interfero- 
metric gravitational wave detectors. 
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IMPLICATIONS FOR ADVANCED DETECTORS 

The low mechanical loss of silica in the 10 - 1000 Hz 
bandwidth, coupled with its optical and thermal prop- 
erties, makes it an attractive material for the optics of 
next generation interferometric gravitational wave detec- 
tors. Fused silica has recently been chosen as the test 
mass substrate for Advanced LIGO Q, which has been 
approved and recommended for funding by the US Na- 
tional Science Foundation. 

If the bulk and surface loss predicted herein can be 
achieved, the mirror thermal noise in Advanced LIGO 
with fused silica mirrors will likely be dominated by the 
coating [6lR. l2ll| . The mirror thermal noise contributions 
to the total Advanced LIGO noise budget are shown in 
FigureEl Table UTTI shows the predicted sensitivity of Ad- 
vanced LIGO with silica optics to two possible sources of 
gravitational waves: binary neutron star inspirals (BNSI) 
and binary 10 Mq black hole inspirals. Two different 
scenarios of coating thermal noise are shown: the best 
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